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ABSTRACT 

We present Submillimeter Array observations of the massive star-forming re- 
gion IRAS 20126-1-4104 in the miUimeter continuum and in several molecular line 
transitions. With the SMA data, we have detected nine molecular transitions, 
including DCN, CH3OH, H2CO, and HC3N molecules, and imaged each molec- 
ular line. From the 1.3 mm continuum emission a compact millimeter source is 
revealed, which is also associated with H2O, OH, and CH3OH masers. Using a 
rotation temperature diagram (RTD), we derive that the rotational temperature 
and the column density of CH3OH are 200 K and 3.7 x 10^^ cm~^, respectively. 
The calculated results and analysis further indicate that a hot core coincides with 
IRAS 20126+4104. The position- velocity diagrams of H2CO (3^ 3-2^ J and HC3N 
(25-24) clearly present Keplerian rotation. Moreover, H2CO (3o3-2q 2) is found 
to trace the disk rotation for the first time. 

Subject headings: ISM: individual (IRAS 20126-1-4104) — ISM: kinematics and 
dynamics — ISM: molecules — stars: formation 
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1. Introduction 

Massive star formation, especially in the earliest phase, is poorly understood, due to 
the large distances involved (> 1 kpc), complex cluster formation environments and shorter 
evolutionary timescales of massive stars. Hot cores represent the early phase in the process 
of massive star-formation (Qin et al. 2010 ). The hot cores are defined as compact ( < 
0.1 pc, n > 10^ cm~^), relatively high temperature (Tk > 100 K) cloud cores (Kurtz et al. 
2000). Thus, observation of hot cores with high angular resolution and various molecular 
lines can provide important information on the early evolutionary phase of massive star 
formation. 

IRAS 20126+4104, at a distance of 1.7 kpc from the Sun (Dame & Thaddeus, 1985; 
Wilking et al. 1989; MoscadeUi et al. 2011), is a massive star-forming region with a 
bolometric luminosity of 1.3xlO^L0. Previous observations of NH3, CH3CN and CS 
molecular lines toward IRAS 20126+4104 indicated that a Keplerian disk was detected 
around a ~ 7-24 Mq point mass (Zhang et al. 1998; Keto & Zhang 2010; Cesaroni et al. 
1999, 2005). Traced in CO, NH3, SiO, HCO+, CH3OH and H2 molecular hues (Shinnaga, 
et al. 2008; Su ct al. 2007; Zhang ct al. 1999; Cesaroni et al. 1997, 1999, 2005; Shepherd et 
al. 2000), a bipolar jet/outflow system was widely investigated toward IRAS 20126+4104. 
In addition, H2O, OH and CH3OH masers in IRAS 20126+4104 were revealed from 
observations at centimeter bands (Tofani et al. 1995; Trinidad et al. 2005; Minier et al. 
2001; Edris et al. 2005). The CH3OH masers might be associated with a hot core (Minier 
et al. 2001). Prom VLA observations toward IRAS 20126+4104, Tofani et al. (1995) did 
not detect the 1.3 and 3.6 cm continuum emission, but Hofner et al. (1999) detected two 
elongated structures at the 3.6 cm continuum band. They concluded that the emission 
was associated with thermal ionized jets and not from a UCHii region. Moreover, IRAS 
20126+4104 has the typical colors of UCHii (Wood & Churchell 1989), suggesting that 
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IRAS 20126+4104 may be a precursor of a UCHii region (Molinari et al. 1996). Thus, 
IRAS 20126+4104 provides us with an opportunity to study the early evolution of massive 
star formation. 

We have carried out multiline observations toward the massive star-forming region 
IRAS 20126+4104 with the Submillimeter Array. Various molecular lines are used to 
investigate the physical and chemical processes of IRAS 20126+4104. So far, only a few 
molecular lines exclusively indicate rotational motions (Beuther 2007). In this paper, we 
find new molecular line ( H2CO ) which can trace the Keplerian disk. In Sect. 2, we 
summarize the observations. In Sect. 3, we give the general results. We discuss core 
properties and core kinematics in Sect. 4.1 and Sect. 4.2, respectively. In Sect. 5, we 
summarize our main conclusions. 



2. OBSERVATIONS AND DATA REDUCTION 

The data are from SMA archive. Observations toward IRAS 20126+4104 were 
carried out with the SMA on 2004 August 23, at 218 (lower sideband) and 228 GHz 
(upper sideband). The typical system temperature was 130 K. The compact configuration 
with seven antennas was used. The phase center was R.A(J2000.0)= 20'^14™26^.00 and 
decl.(J2000.0) = 41°13'31".50. The spectral resolution is 0.8125 MHz, corresponding to a 
velocity resolution of 1.1 km s""^. The absolute fiux density scales were determined from 
observations of Uranus, and the bright quasar 3C 279 was used for bandpass calibration. 
QSO 2015+371 and QSO 2202+422 were observed for the antenna gain corrections. The 
calibration and imaging were performed in Miriad. The continuum map was made using 
the line-free channels. The spectral cubes were constructed using the continuum-subtracted 
spectral channels. Self-calibration was performed to the continuum data. The gain solutions 
from the continuum were applied to the line data. The synthesized beam size of continuum 
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and line images with robust weighting was approximately 1".26 x 0".96 (P.A. = 89°. 0). 

3. RESULTS 

3.1. Continuum Emission at 1.3 mm 

Figure 1 is the 1.3 mm continuum map of IRAS 20126+4104 obtained by combining 
the lower and upper sideband data. The continuum emission shows a compact source, 
which is unresolved in the beam at 1.3 mm. By using a two-dimensional Gaussian fit for 
the continuum emission, we obtained that the total flux density is 0.62 ± 0.03 Jy, the 
deconvolved source size is 0.8" x 0.5" (P.A. = 87.9°), and the peak position is R.A.(J2000) 
= 20^14''^26.'040 (AR.A. = ±0.04"), decl.(J2000) = +41°13'32."55 (Adecl. = ±0.03") 
with an intensity of 0.46 ± 0.05 Jy beam~^. The peak position of the 3.3 mm continuum 
emission is R.A.(J2000) = 20^14"^26.^008 (AR.A. = ±0.01"), decl.(J2000) = ±41°13'32."73 
(Adecl. = ±0.01") (Cesaroni et al. 1997). For the 7 mm continuum emission (Hofner et al. 
1999), the peak position is R.A.(J2000) = 20^14"^26.^026, decl.(J2000) = ±41°13'32."70. 
The three positions coincide within the uncertainties. H2O, OH and CH3OH masers in 
IRAS 20126±4104 have been detected from observations at centimeter wavelengths (Tofani 
et al. 1995; Trinidad et al. 2005; Minier et al. 2001; Edris et al. 2005). In Figure 1, H2O 
masers are shown with the blue flUed triangles. The green open squares present OH masers. 
The flUed black circles show CH3OH masers. Most of the masers are associated with the 
1.3 mm continuum emission, but are not located at the peak position of the 1.3 continuum 
emission. 
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3.2. Molecular Line Emission 

In the frequency coverage of the 2 GHz band, deuterated hydrogen cyanide (DCN), 
formaldehyde (H2CO), cyanoacetylene (HC3N) and methanol (CH3OH) emission lines were 
identified by use of the CDMS and JPL catalog ( Miiller et al. 2005; Pickett et al. 1998), 
as done by other authors in the same frequency range toward Orion (Sutton et al. 1985) 
and Sgr B2 (Qin et al. 2008). The transition, rest frequency and upper level energy (£■„) of 
each molecular line are presented in Table 1. 

Molecular line images can provide valuable information on the spatial distribution of 
molecules and kinematics, then we imaged the molecular lines of DCN, CH3OH, H2CO, and 
HC3N. The 1.3 mm continuum emission has been subtracted before computing these maps. 
Figure 2 presents the integrated intensity maps of each molecular line. " x " indicates 
the peak position of the continuum emission. In Figure 2, each line emission is spatially 
coincident with the continuum emission. 

Spectra of molecular emission are present in Figure 3. The spectra were extracted 
from the peak positions of each integrated intensity map. We have made Gaussian fits 
to all the spectra. The central hne velocity (Vlsr), peak intensity (/p), and full width at 
half-maximum (Ay) are summarized in Table 1. 



4. DISCUSSION 

4.1. Core Properties 

Our interferometric continuum and multiline observations toward the massive star- 
forming region IRAS 20126-1-4104 reveal a compact millimeter source, which is associated 
with H2O, OH, and CH3OH masers. Millimeter continuum emission mainly traces warm 
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dust heated by the embedded young stellar object. Four transitions of CH3OH have been 
detected in IRAS 20126+4104, so we can use a rotation temperature diagram (RTD) to 
estimate the rotation temperature and the column density. Assuming local thermodynamic 
equilibrium (LTE), lines being optically thin and gas emission filling the beam, the rotation 
temperature and the beam-averaged column density can be expressed by (Goldsmith & 
Langer 1999; Liu et al. 2002; Qin et al. 2010) 

H^) = H-^) - ^, (1) 

9u ^rot rot 

Where A'^ is the column density of the upper energy level, is the degeneracy factor in 
the upper energy level, Nt is the total beam-averaged column density, Qrot is the rotational 
partition function, is the upper level energy in K, and T^-ot is the rotation temperature. 
Figure 4 is the rotation temperature diagram. A linear least-squares fit is performed 
toward the four transitions of CH3OH. The RTD can be corrected by multiplying Nul 9u 
by the optical depth correction factor Cr — t/{1 — e'^), where r is the optical depths. 
According to Qin et el. (2010), we iteratively apphed the Cr correction to the RTD until the 
solution converged. After the optical depth is corrected, the rotation temperature and the 
beam-averaged column density are 200 ± 57 K and (3.7 ± 1.3) x 10^^ cm~^, respectively. 
The rotation temperature of 200 ± 57 K is the same as that estimated in the CH3OH (2-1) 
and (5-4) fines (171 K; Cesaroni et al. 2005), and in the CH3CN fines (201 K; Cesaroni et 
al. 1997). 

In the Orion molecular cloud, the strongest CH3OH emission is from the compact 
ridge (Schilke et al. 1997), giving fractional abundances relative to H2 of 1.2 x 10^^. 
We adopt 0.007 pc as the deconvolved size of the CH3OH emission derived from the 
integrated intensity maps of CH3OH, which is smaller than the canonical size of 0.1 pc that 
is commonly quoted for hot molecular cores (Kurtz et al. 2000). By using the fractional 
abundance in the Orion molecular cloud, the calculated H2 density is 4 xlO^ cm~^. The 
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relatively higher gas temperature (200 ± 57 K), H2 density (10^cm~^), and smaller size 
(0.007 pc) indicate a hot core in this region. Following Gibb et al. (2000), we adopt 3 
X 10^^ cm~^ as a typical H2 column density. We obtain that the fractional abundance of 
CH3OH relative to H2 is 8 x 10~^. The fractional abundance of CH3OH is close to that in 
the Orion hot core (Blake et al. 1987), but higher than that in the Sgr B2 core (Qin et al. 
2008). 

4.2. Core Kinematics 

The different characteristics of molecular lines may be used as tools to trace various 
physical processes of massive star formation. To understand the kinematic signatures 
of the rotating disk, we use the H2CO (3^ 3-2^ 2) and HC3N (25-24) lines to construct a 
position-velocity diagram cut through the continuum peak position at a position angle of 
58°. In Figure 5, the continuum emission is redshifted to the southwest and blueshifted 
to the northeast. The position-velocity diagrams of H2CO (3(, 3-2^ 2) and HC3N (25-24) 
clearly present Keplerian rotation. The Keplerian rotation traced in the HC3N (25-24) 
line is extremely similar to the model disk constructed by Zhang et al. (1998). So far, for 
many sources only a few molecular lines exclusively indicated rotating motions in massive 
star formation (Beuther 2007), e.g., CH3CN, C^^S, NH3, HCOOCH3, C^^O, H^^O, HN13C, 
HC3N and CH3OH. Here we find a new molecule (H2CO 3^ 3-2,32), which can be used to 
investigate the rotational signatures of a disk. 

In Figure 5, a velocity gradient in the NE-SW direction is clearly depicted. The 
velocities of the two emission peaks are -6.0 and -2.0 km s~^ with 0.60" separation for 
the H2CO (3q 3-2q 2) molecule, and -5.6 and -3.0 km s~^ with 0.55" separation for the 
HC3N (25-24) molecule, respectively. Assuming equilibrium between the rotational and 
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gravitational force at the outer radius of the Keplerian disk: 

Mrot = ^ = 2.3 X 10^ V\\im s-i)r(pc) Mq, (2) 

where V is the velocity difference, r is the spatial separation of the emission peaks and G 
is the gravitational constant. The derived dynamical mass is 18.0 Mq and 7.1 Mq in the 
H2CO (3(j 3-2(, 2) and HC3N (25-24) molecular hues, respectively. The calculated dynamical 
mass is smaller than a mass of ~ 20 Mq obtained by Zhang et al (1998) in the NH3 
molecular lines. 

5. SUMMARY 

We have observed the massive star-forming region IRAS 20126-1-4104 in the millimeter 
continuum and in several molecular line transitions. The 1.3 mm continuum emission 
reveals a compact millimeter source which is associated with H2O, OH, and CH3OH masers. 
From the SMA data, nine molecular transitions, including DCN, CH3OH, H2CO, and HC3N 
molecules are detected. The rotational temperature and the column density of CH3OH 
emission are 200 K and 3.7 x 10^^ cm"^, respectively. The calculated results and analysis 
further indicate that a hot core coincides with IRAS 20126+4104. A Keplerian rotational 
disk in the IRAS 20126+4104 hot core is verified in the H2CO {\^-\^) and HC3N (25-24) 
lines. H2CO (3(, 3-20 2) is found to trace the disk rotation for the first time. The dynamical 
mass of the Keplerian rotational disk is 18.0 Mq and 7.1 Mq, determined in the H2CO 
(3(, 3-2(, 2) and HC3N (25-24) molecular lines, respectively. 

We thank the SMA staff for the observations and the anonymous referee whose insights 
greatly improved this manuscript. We also thank Dr. H. Shi and other people for providing 
help with the data reduction and analysis. This work was supported by the National 
Natural Science Foundation of China under Grant No. 10473014. 
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20''14"26f3 26=2 26?1 26?0 25?9 25?B 
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Fig. 1. — Continuum map toward IRAS 20126+4104 at 1.3mm obtained by combining the 
data from both sidebands (218 and 228 GHz). The contours are at -3, 3, 4, 5, 6, 7, 8, and 9 a. 
The rms noise level is 0.05 Jy beam^^(lcr). The synthesized beam is 1.3" x 0.8", P.A.= 88.9° 
{lower right corner). " x " indicates the peak position of the continuum from the Gaussian 
fit. H2O masers are shown with blue filled triangles. The green open squares present OH 
masers. The filled black circles show CH3OH masers. 
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Fig. 2. — (a) The DCN 3-2 integrated intensity map. The contour levels are at -2.2, 2.2, 4.4, 
6.6, 8.8, 11.0, 13.2 and 14.3 Jy beam-^ km s'^. (b) The H2CO 3o,3-'^o,2 integrated intensity 
map. The contour levels are at -1.0, 2.5, 5.0, 7.5, 9.0, 12.5, 15.0 and 16.2 Jy beam~^ km s~^. 
(c) The H2CO 32 1-22 integrated intensity map. The contour levels are at -1.5, 2.4, 4.8, 7.2, 
9.6, 12.0, 14.4, 16.8 and 18.1 Jy beam"^ km s~^ (d) The HC3N 24-23 integrated intensity 
map. The contour levels are at -1.5, 2.3, 4.5, 6.7, 8.9, 11.1, 13.3 and 15.1 Jy beam~^ km 
s~^. (e) The HC3N 25-24 integrated intensity map. The contour levels are at -1.6, 2.1, 4.2, 
6.3, 8.4, 10.5, 12.6, 14.7 and 15.9 Jy beam-^ km s-^ (f) The CH3OH 5^^-422 integrated 
intensity map. The contour levels are at -1.1, 2.2, 4.4, 6.6, 8.8, 11.0, 13.2, 15.4, 17.6, and 
19.3 Jy beam"^ km s~^. (g) The CH3OH 5-^^-4^^ integrated intensity map. The contour 
levels are at -2.1, 2.8, 5.6, 8.4, 11.2, 14.0, 16.8, 19.6 and 22.8 Jy beam"! km s'^. (h) The 
CH3OH 21i 2o-21o,2i integrated intensity map. The contour levels are at -1.4, 1.8, 3.6, 5.4, 
7.2, and 8.2 Jy beam~^ km s~^. (i) The CH3OH 16^15-15213 integrated intensity map. The 
contour levels are at -1.5, 2.0, 4.1, 6.1, 8.2, 10.2, 12.2, 14.3 and 16.1 Jy beam"! km s'^ " x " 
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Table 1: Observed Parameters of Each Line 



Molecule 


Transition 


Rest Frequency 


En 






AV 


Channel rms 






(GHz) 


(K) 


(km s"-*^) 


(Jy beam"^) 


(km s""*^) 


(Jy beam~^) 


DCN 


3-2 


217.239 


20.9 


-3.4(0.2) 


1.65(0.08) 


8.0(0.4) 


0.2 


H2CO 


3 -9 

"0,3 0,2 


218.222 


21.0 


-3.8(0.3) 


2.38(0.19) 


8.6(0.8) 


0.3 






218.760 


68.1 


-5.1(0.3) 


2.09(0.13) 


7.5(0.6) 


0.2 


HC3N 


24-23 


218.325 


131.0 


-4.1(0.2) 


1.78(0.11) 


7.8(0.6) 


0.2 




25-24 


227.419 


141.9 


-4.6(0.1) 


2.53(0.09) 


6.7(0.3) 


0.2 


CH3OH 


5 -4 


216.946 


55.9 


-3.8(0.2) 


2.02(0.12) 


8.2(0.5) 


0.2 




4 -3 

^2,2 '-'1,2 


218.440 


45.5 


-4.1(0.1) 


3.01(0.17) 


7.3(0.5) 


0.2 




21 -21 

"^-•-1,20 ^-"-0,21 


227.095 


557.1 


-3.7(0.4) 


0.87(0.06) 


9.3(0.9) 


0.1 






227.815 


327.2 


-4.0(0.2) 


1.62(0.07) 


8.2(0.4) 


0.1 
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Fig. 3. — Molecular spectra at the peak position of different molecular core. The solid and 
dashed curves are the observed spectra and the Gaussian fitting to the spectra, respectively, 
(a) DON 3-2. (b) H2CO \^-2^^^. (c) H2CO \^-2^^^. (d) HC3N 24-23. (e) HC3N 25-24. (f) 
CH3OH 5,,,-4,^. (g) CH3OH 4,_,-3,,,. (h) CH3OH 21,^„-21„_,,. (i) CH3OH 16,_,,-15,,,3. 
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Fig. 4. — Population temperature diagram of the observed CH3OH transitions. The vertical 
bars mark the la(N^/g^) errors from the integrated intensities. The linear least-squares fit 
(solid line) gives a rotation temperature of 200 ± 57 K. 
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(a) HjCO 3o.3-2o.2 




(b) HCjN 25-24 
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Fig. 5. — (a) The position- velocity plot of the H2CO 3-2;, ^ emission at a position angle of 
58° along the major axis of the disk (Zhang et al. 1999). The contour levels in observed data 
are plotted at every 2a. la = 0.6 Jy. Position offsets is (0.6", 1.0"). (b) The position-velocity 
plot of the HC3N 24-23 emission at a position angle of 58° along the major axis of the disk. 
The contour levels in observed data are plotted at every 2a. la — 0.5 Jy. Position offsets is 
(0.6", 1.2"). The horizontal dashed hne marks the center of the disk, and the vertical dashed 
line marks the cloud systemic velocity. The full lines show Keplerian rotation curve. 



